
subroutine compute_Richter_corr(dst,corr,nv,RichDist,RichVal)
! ATTENZIONE: questa routine funziona solo se la tabella
! inizia dalla distanza 0.

integer dst
real*4  corr,dd
integer nv,i,RichDist(100)
real*4       RichVal(100)

dd=float(dst)

do i=1,nv

  if(RichDist(i) .gt. dd)then
    a = RichDist(i) - RichDist(i-1)
    b = RichVal(i) - RichVal(i-1)
    corr = RichVal(i-1) + b/a*(dd - RichDist(i-1))
    return
  endif 
enddo
! oltre la distanza massima:
corr = RichVal(nv)
return
end
    